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ABSTRACT 

An analytic model is proposed for non-radiating accretion flows accompanied 
by up or down winds in a global magnetic field. Physical quantities in this 
model solution are written in variable-separated forms, and their radial parts are 
simple power law functions including one parameter for wind strength. Several, 
mathematically equivalent but physically different expressions of the criterion 
for wind generation are obtained. It is suggested also that the generation of 
wind is a consequence of the intervention of some mechanism that redistributes 
the locally available gravitational energy, and that the Bernoulli sum can be a 
good indicator of the existence of such mechanisms. 

Subject headings: accretion, accretion disks — magnetohydrodynamics: MHD 
— galaxies: nuclei, jets 
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1. INTRODUCTION 

Blandford & Payne (1982) have derived a criterion for the generation of centrifugal 
winds from magnetized accretion disks. It has been shown in their self-similar solution that 
a centrifugal wind appears if the inclination of the poloidal magnetic field lines penetrating 
a Keplerian accretion disk makes an angle of less than 60° with the disk surface. Although 
this result is intuitively very understandable, in a more general situation that does not 
satisfy self similarity, this kind of criterion becomes meaningless because the inclination 
angle may be different at different locations even along a stream line. Therefore, we need 
to seek for physically more essential expressions of the criterion for the generation of winds 
from accretion disks, in order to improve the understanding of such processes (for general 
reviews of the wind and jet theories, see e.g., Begelman, Blandford & Rees 1984; Ferrari 
1998; Livio 1999). 

In relation to recent development of the theories of optically thin ADAFs (advection- 
dominated accretion flows: see e.g., Ichimaru 1977; Narayan & Yi 1994, 1995; Abramowicz 
et al. 1995; Chen et al. 1995), the sign of the Bernoulli sum has drawn considerable 
attention (and has caused also confusion) as a possible indicator of the presence of 
unbounded outflows such as winds or jets. The sum consists of the gravitational energy, 
kinetic energy and enthalpy, all per unit mass of a fluid. 

Narayan & Yi (1994, 1995) have shown that the sign of the Bernoulli sum is necessarily 
positive in their self-similar ADAF solution, and argued that this is a genuine property 
of general ADAFs. They have suggested also that a flow with positive BernouUi sum 
(and hence ADAFs, specifically) can easily drive winds and jets. Further developed and 
popularized by Blandford & Begelman (1999) in their influential "ADIOS" (adiabatic 
inflow-outflow solution) paper, these suggestions have been widely accepted. The latter 
paper stresses and demonstrates the necessity for including winds or jets in constructing a 
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satisfactory theoretical model of inefficiently radiating accretion flows. 

Meanwhile, several authors have demonstrated that the positivity of the Bernoulli sum 
is not a genuine property of ADAFs. As shown by Nakamura (1998), the analytic solution 
of Honma (1996) that includes a diffusion cooling by turbulence has negative sign at least 
in a special case of 7 = 5/3, where 7 is the index of polytrope. Further it has been shown 
that, under the influence of convective coohng (Igumenshchev & Abramowicz 2000) and of 
the inner and outer boundary conditions (Abramowicz, Lasota, & Igumenshchev 2000), low 
viscosity ADAFs (with a < 0.3, where a is the viscosity parameter) have negative sums. 
This result agrees with all of the subsequent 1-D, 2-D and 3-D numerical simulations by 
Abramowicz and his collaborators. For inviscid ADAFs, the no-wind solution in a global 
magnetic field (Kaburaki 2000, hereafter KOO) shows that the Bernoulli sum is always zero 
within the adopted approximation. 

Another point to be mentioned is the claim by Abramowicz, Lasota, & Igumenshchev 
(2000) that a positive Bernoulli sum is only a necessary but not the sufficient condition for 
unbounded outflows, as shown for inviscid, non-magnetic fluids. They also argue that the 
best example in which a positive Bernoulli sum does not imply unbounded outflows is the 
classical Bondi solution for spherical accretion flows. 

To summarize the present status of the Bernoulli sum described above, at least, the 
positivity of the sum does not seem to be a genuine property of ADAFs. The sign may 
depend on various conditions such as dissipation mechanisms (viscous or resistive), the 
presence of energy transport in a fluid (convection, conduction etc.) and the effects of 
boundaries. However, it is still controversial whether the sign of the Bernoulli sum can be a 
direct indicator of the presence of unbounded outflows. 

It seems natural to believe that the sum is everywhere zero (the asymptotic value 
at inflnity) in a fluid, as far as the radiation cooling and other mechanisms of energy 
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redistribution in the fluid are completely negligible, since then all the dissipated energy 
remains as thermal energy of each fluid element and is merely advected down the flow (i.e., 
the flow is completely advective). Therefore, a non-zero Bernoulli sum should indicate the 
intervention of some mechanism of energy redistribution within a non-radiative flow. We 
may expect convection, conduction or fluid viscosity as such mechanisms of redistribution. 

In the present paper, we will derive an analytic model that describes non-radiating 
accretion flows accompanied by up or down winds in a global magnetic fleld, based on our 
previous treatment, KOO. Although the main motivation of the present work is the desire 
to improve broadband spectral fittings to low-luminosity active galactic nuclei (LLAGNs) 
and normal galaxies (e.g. for Sgr A*, see Kino, Kaburaki & Yamazaki 2000; Yamazaki, 
Kaburaki & Kino 2001), the resulting solution will be useful to discuss the issues about the 
Bernoulli sum discussed above. 

2. SIMPLIFYING ASSUMPTIONS 

A schematic drawing of the global configuration presumed in the present consideration 
is given in Figure 1 of KOO. The viscosity of accreting plasma is completely neglected to 
clarify the role of magnetic field. An asymptotically uniform magnetic field is vertically 
penetrating the accretion disk and is twisted by the rotational motion of the plasma. 
Owing to the Maxwell stress of this twisted magnetic field, a certain fraction of the angular 
momentum of accreting plasma is carried out to infinity, and this fact enable the plasma 
gradually infall toward the central black hole. 

As shown in Appendix of KOO, the geometry of accretion fiows is very essential in 
specifying their physical properties. For example, the ADAF solution follows almost 
automatically from the assumption of constant opening angle of the disk (i.e., A — const.. 



-6- 



where A is the half-opening angle), which respects the spherical nature of the gravity. 
In this sense, spherical polar coordinates (r, 6, (p) are convenient for the discussion of 
ADAFs. For the relevant physical quantities, we follow the notation of KOO unless specified 
otherwise. 

Further simplifying assumptions made in KOO in obtaining the no-wind ADAF solution 
from the set of resistive MHD equations were those of, i) stationarity {d/dt — 0), ii) 
axisymmetry {d/d(p — 0), iii) geometrically thin disk, iv) weakly resistive disk, v) dominance 
of midplane, and vi) no wind. Among these assumptions, only the last one will be removed 
in the present paper in order to obtain the ADAF solution including winds. 

The third assumption means that A <^ 1, and is most effective to simplify the basic 
equations. It has been shown in KOO that the presence of an external magnetic field Bq 
guarantees the reahzation of such a thin disk even in a hot ADAF situation. Reflecting 
this localized structure, we introduce an angular variable rj = {9 — -k/2)/ A. Then it 
becomes clear that a differentiation with respect to 9 gives rise to a quantity of order A~^ 
{d/d9 — A~^d/drj). We can also safely approximate in the disk as sin^ ~ 1 and cos^ ~ 0. 

The statement that a disk is weakly resistive implies that the "characteristic" magnetic 
Reynolds number 3?, whose deflnition will be introduced later, is large in the sense that 
3?^(r) » 1 in the disk except near its inner edge Tin where 3?(rin) = 1. Actually, the terms 
of 0{^~^) were neglected in KOO. In these situations, the externally given magnetic field 
Bq is largely stretched by the rotational and infalling motion of the accreting plasma so 
that the deformation in the disk becomes much larger than the seed field (if we divide 
the total magnetic field in the form B = Bq + b, then |b| ^ |Bo|). The definition of the 
disk's outer edge is the radius within which the deformation of the poloidal magnetic field 
becomes noticeable (i.e., |bp| ~ |Bo|). 

The assumption v) naturally follows from the consideration that, since majority of 
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matter is concentrated around the midplane of the accretion disk, its physical properties 
should be controlled mainly by this part of the disk. Therefore, we may seek the 
approximate solution that is accurate near the midplane even if making a sacrifice of the 
accuracy at its upper and lower surfaces. According to this spirit, we ignore the quantities 
that are proportional to tanh^Ty since tanh^Ty <^ 1, sech^jy. 

The last condition in the above list is, in fact, not indispensable in obtaining a resistive 
ADAF solution in a global magnetic field (or resistive ADAF solution). Indeed, an accretion 
flow accompanied by a converging flow (or down wind) toward the disk surface has been 
obtained by the present author (Kaburaki 1987) omitting this assumption. It was imposed 
again in KOO, however, in order to understand the energy budget clearly and to firmly 
estabhsh a basic analytic model in the resistive ADAF regime. This condition actually 
consists of two equations, = 0, and jo = for consistency. The former results in an 
r-independent mass accretion rate M, and the latter specifies the radial dependence of the 
toroidal magnetic field as 6^ oc 1/r since je — — (c/47rr) d{rbip)/dr under the assumption of 
ii). 

In relation to the no- wind resistive ADAF solution obtained in KOO, it should be 
stressed that the result actually contains the effects of finite resistivity to the first order 
in the smallness parameter 3fJ~^. The ratio of poloidal to toroidal magnetic field is small 
(i.e., br/bip ~ 3?~^) reflecting the strong twisting of the poloidal seed field by the rotational 
motion. This ratio is maintained by the balance between this twisting and untwisting by 
resistive diffusion. The thin disk structure of the fiow is maintained by the vertical force 
balance between the magnetic pressure of toward the equatorial plane and the opposing 
gas pressure in the disk. Reflecting this plasma enhancement in the disk, the gas pressure 
and density in the disk are the quantities of order 3?^. The finite thickness of the disk 
itself is also a consequence of non-zero resistivity (indeed, A oc 3?"^). The infall velocity is 
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small in the sense that v^jv^p ~ 3?"^. Since the magnetic Reynolds number is a function 
of r, the solution does not have self similarity (i.e., in addition to the presence of angular 
dependences, both the ratios h^jh^p and Vrjv^ vary with the radial distance). 



In order to obtain a more general form of the resistive ADAF solution, we remove here 
the no-wind condition from our list of simplifying assumptions. Fortunately, however, it 
turns out that all but one terms can be omitted finally among the newly appeared terms in 
the leading order equations in A, by the aid of the assumptions iv) and v). We shall check 
this point below, expecting that the angular dependence of every quantity (see KOO, and §4 
below) is not affected by the removal of no-wind condition. 

The component expressions of the resistive MHD equations simplified under the 
assumptions i) through iii) are as follows. They are correct to the leading order in the 
powers of A. In deriving these equations, all quantities except 6e and vq (which are of the 
order of A as confirmed by equations @] and 0) are regarded as of order unity with respect 



3. 



REMOVAL OF NO- WIND CONDITION 
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• equation of motion 
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(y9-component 



dr Ar dr] I r 



(3) 



induction equation 

poloidal component 



1/7 1 \ c 1 dbr 

- {vrbe - vebr) = — - — (4) 

c AnaA r or] 



y3-component 

r%^l. (^^) _ l(ru,b ) + ^ = (5) 

^ dr \ r J dr '" AnaA'^r dri"^ 



mass continuity 



1 (9 , 9 , 1 (9 , , 

- jr{r^P^r) + -^ j^{pve) = Q (6) 
r^ or Ar or] 



• magnetic flux conservation 

In r-component of the equation of motion, there appear two new terms reflecting the 
removal of the no wind condition. They are the second term on the left-hand side and the 
last term on the right. Both of them are dropped, however, from this equation because of 
the assumption v) since they are proportional to tanh^ r]. Further, the first term on the 
left and the third term on the right are dropped, as in the no-wind case, owing to the 
assumption iv) since Vr oc 3?^^ and p oc 3?^^. Thus, the equation becomes finally so simple 
as 

GM v% 1 dv 

— 2~ = TT- y^i 

J. p Or 

There is no new term in equation (Q), and the second term on the left is neglected 
owing to the assumption iv). The resulting equation describes the magnetic confinement 
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of the disk plasma by b^. In (y9-component of the equation of motion, the second term on 
the left and the first term on the right contain vq and jg, respectively. However, the former 
vanishes since f ^ is independent of 6 and the latter drops because of tanh^ rj dependence. 
The resulting form is, as before, 

^^a(riv) ^ be dkp 
dr AirpA dt] 

The poloidal components of the magnetic induction equation are degenerate reflecting 
the degeneracy in Maxwell's equations. The second term on the left of this equation can be 
dropped again by the assumption v), and we have 

1 db 

There appears no change in equation (|5|). 

Thus, the only change that cannot be dropped in the set of basic equations describing 
magnetized ADAFs including winds is the second term on the left-hand side of equation 
(^), which results in a radius-dependent mass accretion rate. The equation of magnetic flux 
conservation is the same as before. 



4. SEPARATION OF VARIABLES 

As written out in KOO, the the set of resistive MHD equations ([1] ~ [4] there) are 8 
equations for 8 unknowns (i.e., for p, p, v and B), and hence the set is closed apparently. 
Other quantities such as j, E, q (the charge density), and T (the temperature) are calculated 
from the subsidiary equations, (5) and (6) in KOO. Actually, however, the above main 
set is not closed because of the degeneracy in Maxwell's equations. Usually, this point is 
supplemented by adding one more equation relating the transfer of energy. 

Although, in many cases, the polytropic relation is adopted for simplicity as such 
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an equation, its validity is rather doubtful except in the special cases of adiabatic and 
isothermal processes. On the other hand, the full equation for the energy transfer, not only 
in the fluid but also to the surroundings by the radiation, is too complicated to be treated 
analytically. Here, we would rather let the set being open as in KOO since an open set does 
not mean the absence of solutions. Actually, it means only that there is no systematic way 
of solving it. 

Instead of adding energy equation, we have put two specific constraints in KOO in order 
to obtain the ADAF solution without wind. Before entering into the discussion of them, we 
first separate the variables in the following form being led by the experience in KOO. 



br{^,v) = br{^) sech^Tjianhrj, (11) 

be{C,v)= 6e(e) sech^r;, (12) 

Ki^,v) = -KiO^^'^^v, (13) 

Vr{C,v) = -VriO sech^77, (14) 

f^e(^,^) = ^e(Otanh77, (15) 

v^itv) = M0^ (16) 

P{C,V) = P(0 sech^r?, (17) 

p{^:V)= P(0sech\ (18) 

TiC,v)-m, (19) 

Jr-(e,^) = -Jr(0sech2,y, (20) 

je{i,v)=MOts^iahr], (21) 

J^(e,r?) = -i^(Osech\ (22) 

£;,(e,77) = £^.(e)sech\ (23) 
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EJ^,V) = EM) sech^77, 



(24) 
(25) 



where the radial coordinate is normahzed by a reference radius ro as ^ = r/rg. In the 
problems of accretion in a asymptotically uniform magnetic field, it is natural to choose rout 
as To- The sign of veiCv) is chosen for positive vg to correspond to an up wind (outflow) 
from the disk surfaces. 

The set of basic equations are now rewritten as the set of ordinary differential equations 
for the radial part functions: 



r 



+ 



1 dp , GM 
p dr 

P 



dl 



Vr— — 



Stt' 
1 heh^ 



dr 

Vrbe 



AttA P 



47r(7A r ' 



,dn 
dr 



^ = 0, 



Vr 



A 



r-^ In^r'^pVr] 
dr 

r— ln(r%r) 
dr 



(26) 

(27) 
(28) 
(29) 
(30) 
(31) 
(32) 



where we have defined I = rv^ and fl = v^/r. The angular dependences of the relevant 
quantities are mutually consistent within the assumption v). 

In addition to the above set, the subsidiary equations are 



T 



fj_ p 
R p 



47rA r ' 



(33) 
(34) 
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^. = 4^(r^.), (37) 

-Eg = ^ {Vrh^ - vj)r) , (38) 

where R is the gas constant, p, is the mean molecular weight, and should vanish owing 
to the assumption ii). 



5. IRAF CONDITIONS 

As mentioned in the previous section, two specific constraints have been placed in 
obtaining the solution in KOO. Then, the solution has been retrospectively shown to be 
fully advective. Thus, the constraints have taken the place of energy equation, and selected 
the ADAF solution among others. Therefore, they may be called the ADAF conditions. 

We shall adopt the same constraints also in the present paper in order to select the 
corresponding specific kind of solution to the above set of resistive MHD equations. In 
view of the resulting solution, however, it seems to be needed to extend the concept of 
ADAF here. The essence of the conventional ADAF is in that the flow is a very ineffective 
emitter of radiation field, rather than in the point that internal energy is transported 
mainly in the form of advection. Such a character is therefore better specified by the term, 
inefficiently-radiating accretion flow or "IRAF". In such a flow, thermal energy may in 
general be transported by convection or by conduction as well as by advection. Therefore, 
we call hereafter the specific constraints the IRAF conditions. 
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The first of tliem is 



a 



1 dp I GM 



pdr I 
1 dp I GM 



const. (< 1). 



(40) 



p dr / 



Tlien, equation (P) or ( 126|) results in a reduced Keplerian rotation of the form 



(41) 



where v\<^{r) = (GM/rY^"^ is the Kepler velocity. The first IRAF condition requires that the 
pressure gradient term, which appears in the equation of radial force balance (|) together 
with the centrifugal force density, behaves exactly in the same way as the gravitational force 
density (i.e., they have the same r- and ^-dependences). In general, however, the former 
is determined as a complicated consequence of the thermal processes such as dissipative 
heating, radiative cooling and advective cooling, and hence may have different r- and 
^-dependences from the latter. In spite of this affair, equation (^) demands the pressure 
gradient to adjust itself to follow the form of the gravity. Here, one can see the dominance 
of the gravity over the thermal processes in the IRAF solution. 

The second IRAF condition requires that 



Similarly to the first condition, the gravity that specifies the rotation law in the denominator 
determines the degree of magnetic field twist as expressed by the numerator. In general, 
however, the latter is determined by a complicated consequence of the electro dynamic 
processes and hence may have different r-, and ^-dependences from the former. Here, we 
see the dominance of the gravity again, this time, over the electrodynamic processes. 




(42) 
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Then, equation (|), or (|30D reduces to 



47r(l - /3)aA2 
Substituting equation ( 14T|) for Q, we obtain 



(ir 



1 -1 



(43) 



3.(1 -/^)AV 



6^ c 



(44) 



where /k = V GMr. This ratio, which is specified by the reduced-Keplerian angular- 
momentum distribution, is adopted as the characteristic magnetic-Reynolds number of this 
solution. It must be noted in this regard that the actual magnetic Reynolds number that is 
defined by the ratio of the convection to diffusion terms in the magnetic induction equation 
is everywhere and always 1 since, in a stationary state, they are balanced exactly by each 
other. 



6. SOLUTION WITH WINDS 

Although a fully advective solution has been automatically obtained in KOO starting 
from the toroidal field of the form b^p oc which is the consequence of je = 0, we cannot 
use this relation here in obtaining a more general situation with non-zero vg and jg. Instead, 
we start from the following form for the poloidal magnetic field, 

briO = Bo r^'/'-"\ (45) 

where the parameter n specifies the strength of a wind with n = corresponding to the no 
wind case. 

Substituting the above expression into equations ( ^if ) and (|32|) , we obtain 

= 3}„B„ (46) 
«o ^ ^^(^--^'J-ff^"^ (47) 
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and 



beiO = '^^^Bor^'/'-\ (48) 



respectively. Then, it follows from equation p7| ) that 



Pio = ^ r^(^-"^. (49) 

The subscript is referred to each quantity at tq. 

Combining the expressions for br and bg, we obtain from equation (^^) 

^ 3(1 - ami -0) ^ , 

2n + l % 

and further from this result and equation (|3T|), 

in + 1 Jfto 



The density is calculated from equation (|2^) as 

^^^^"24vr(l -«)(!-/?) <o ^ ^ ^ 



With the above expressions for various quantities, we can determine first (3 from the 
second IRAF condition ( ^2]) and then a from the first condition ( ^0]) as 



« = ^(l-n), /5 = ^(l-n). (53) 
Therefore, the final forms of the above tentative expressions are 

m)-%e'\ ^, = {?^fi^ (54) 



-.(0 = y^^e , (55) 
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5.(0 = \/^-Kor^^ (57) 

p(0 = r(^-^"\ (58) 



and the subsidiary equations yield 



KO 



no = I % r\ (59) 



C ^ 3?o5o ^-(2-n) 



C \ 3fJoA) £.-(2-n) 



= = -B^(«) = 0. (63) 

It turns out from the above expressions that the resistive IRAF solution has some 
characteristic features. First of all, temperature is independent of n and always has the 
virial value. The radial powers of the all components of velocity are independent of n, 
while their coefficients depend on n. On the other hand, the powers of other quantities 
are dependent on n, while their coefficients are independent of n except for be and je. The 
vanishing of the electric field reflects one of our implicit assumption that the electric load 
such as the acceleration of a bipolar jet is completely neglected in obtaining the solution, 
and hence the total electric power available is spent in the accretion disk. 



7. ENERGY BUDGET 

Some global considerations of the physics associated with our accretion disk reveal 
more information about the present solution. First, from the deflnition of the disk's inner 
edge (i.e., 3fJ(^in) = 1) and equation (^H), we obtain 

^in = (64) 
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and the condition <^ rout guarantees <^1. 

The total magnetic flux penetrating the disk surface can be regarded as being generated 
by the toroidal current in the disk, since the contribution of the seed field is negligible there. 
The whole flux closes beyond the outer edge, reducing the magnetic flux in that region. 
Then, the reduced amount of flux should appear in the narrow central region within the 
inner edge, with an enhanced strength. Approximating the mean intensity of the central 
magnetic field by 6r(6n), we have 

vrrij, 6r(nn) = / be 2Txr dr, (65) 

which results in the relation 

^n+l/2 

A = ^ C^'/' = (66) 

^ Sin 

The approximate expression holds because ^^^^^'^ ^ 1 (this point will be confirmed at the 
end of this section). 

Owing to the presence of vertical flows, the mass accretion rate becomes radius 
dependent like 

M(0 = - / 27rpyA df] = Mo (67) 



oo 



where 



M„^,^&M1. (68) 



3 VGM 



Solving equation (|^) for tq, we obtain 



It is quite easy to evaluate the Bernoulli sum for the present solution: 

M^,V) ^ l^J-^ + h = ^ v^,r' = Be(0, (70) 
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where h denotes the specific enthalpy and h = {5/2)p/p for any ideal gas. The poloidal 
velocity does not appear in the kinetic energy term, since its contribution is small compared 
with that from the toroidal component by 3fJo ^. The sum is not a constant along each 
stream line because the flow is dissipative. For non-zero n, its absolute value increases 
from zero at infinity as the radius decreases. It is apparent that, within the framework 
of the present solution, the sign of the Bernoulli sum is indeed a discriminator of wind 
generation since it directly reflects the sign of n: i.e., according to whether it is positive, 
zero, or negative we have up wind, no wind, or down wind, respectively. As discussed in §1, 
the appearance of non-zero sum may indicate the intervention of some other mechanism of 
energy transport than advection. This point becomes clear soon. 

The local energy budget is as follows. The Joule heating rate is calculated as 

i2 -2 

a (7 



2n 



1 ^2n+l GMMo ^-2{2-n)gg^i^4^^ (7^^ 



167r 4 

where c^/crA^ and have been eliminated by the aid of equations (15^) and (|68D. The 
advection cooling, on the other hand, is 

Cdv(^' V) = div(/ipv) - (v • grad)p 

— — (r^hpv ) + — — (hpve) - v —- 
An 



+ 1 5j2n+l GMM^ ^^^^4^^ (72) 



167r 

where the last term in the second line should be dropped on account of the assumption v). 

It is evident, therefore, that the advective cooling is in general not balanced by the 
Joule heating alone. This means that there must be some additional heating that have not 
been mentioned explicitly, so that the energy balance in a stationary state should in fact be 

(it + 9add = ?adv (73) 



- 20 - 



since radiation cooling is negligible {q^,^^ ^ ^i-ad) ^^^o in the present solution. From the 
above expressions for gj" and q~^^, we have 



This is what has been suggested just above by the appearance of a non-zero Bernoulli sum. 

Plausible examples of such additional heating mechanisms may be viscous heating (g^s; 
as is popular in the standard and viscous ADAF theories), convection heating {q'^onv, e.g., 
Narayan, Igumenshchev, & Abramowicz 2000; Quataert & Gruzinov, 2000), conduction 
heating (qcond) Among these examples, the first one should always be positive, 

while the others may be negative as well. The above result for q^^^ means that we have up 
wind, no wind or down wind according to whether the additional heating rate is positive, 
zero or negative, respectively. In order to keep the global energy budget in a disk, the role 
of these additional mechanisms should be the redistribution of locally available gravitational 
energy. 

Finally, we derive the restriction on the rage of the wind parameter n. Upper limit 
for n is derived from the radial dependence of p in equation (pH]). Since the density should 
decrease outward, we have n < 1/2. The most stringent lower limit for n is obtained if we 
require that the advective cooling should be really cooling (i.e., — 1/4 < when q~^^ > 0). 
This is equivalent to the requirement that entropy should decrease outward, because in 
a stationary state we can show the equality q~^^ = pT{y ■ grad)s, where s is the specific 
entropy. Thus, we have obtained 




r'^'~"^sech^r/. 



(74) 



Q'add 



1 1 



(75) 




This inequality guarantees 



.n+l/2 



^ 1 as far as ^ 1. 



m 
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8. SUMMARY AND CONCLUSION 

We have extended our former analytic solution for ADAF in a global magnetic 
field to include winds from the disk surface. The solution contains one parameter n 
(—1/4 < n < 1/2) that specifies the sense and strength of the winds. According to whether 
n is positive, zero or negative, we have up wind, no wind or down wind, respectively, with 
the wind strength increased according to its absolute value. 

The following physically different statements are all equivalent mathematically, and 
represent the criterion for the generation of up winds (i.e., n > 0). 

(1) Poloidal magnetic field in the disk decreases with radius slower than r~^^^. 

(2) Rotational velocity of the disk is sufficiently large: i.e., v^p > (1/-\/3)vk- 

(3) Pressure gradient is sufficiently small compare with gravity: i.e., a < 2/3. 

(4) The BernouUi sum in the disk is positive. 

(5) There are other mechanisms of heating that has not been treated explicitly in the 
present consideration. 

Judging from the radial dependence of the poloidal magnetic field, the converging flow 
of Kaburaki (1978) corresponds to the extreme case of n = — 1/4 (down wind). Although 
the poloidal field of the Blandford-Payne solution (1982) formally corresponds to the case 
of n = 1/4 (up wind), the correspondence should not be taken so seriously because their 
self-similar solution is somewhat artificial and does not belong to the same class as the 
solution derived here. The second statement above suggests that the winds are centrifugally 
driven, and the third, that strong winds apt to appear in high pressure environments as 
discussed by Fabian & Rees (1995). 

The fourth and fifth statements are closely related. The BernouUi sum works as a 
discriminator of wind generation, at least within the framework of present treatment. The 
appearance of non-zero sum seems to indicate the intervention of some other mechanisms of 
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energy transport than advection. These additional processes act to redistribute the locally 
available gravitational energy, within a disk. It may be heat conduction or convection that 
carries thermal energy from inner regions to outer . Another important possibility may be 
the viscosity that has been completely neglected in this paper for simplicity, but is generally 
known to have such a redistribution effect (e.g.. Prank, King, & Raine 1992). In either case, 
a more satisfactory solution should be obtained by including each process explicitly in the 
treatment, and it is of course beyond the scope of this paper. 

The wind velocity described by our solution is sub-critical even if it originates from 
an inner portion of the accretion disk, since vq/vk ~ A/3fJ(r) (A -C 1 and K~^(r) < 1), 
and the sound and Alfven velocities are always of the order of the Kepler velocity, i.e., 
Vs(r) ^ Vx(r-) = \j2/3 VKir) where = dp/dp and Vl = 6^/47rp. Therefore, the simple 
coUimation of this type of wind from the surface of an accretion disk cannot explain 
the formation of a relativistic jet. The generation cite of the latter should probably be 
attributed to the central region within the inner edge of disk. In spite of low velocities 
expected for winds from the disk surfaces, their effects on the accretion process can be 
appreciable as seen from the resulting radius-dependent mass accretion rates. 
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